Distillation of mixed-state continuous- variable entanglement by photon subtraction 
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We present a detailed theoretical analysis for the distillation of one copy of a mixed two-mode 
continuous- variable entangled state using beamsplitters and coherent photon-detection techniques, 
including conventional on-off detectors and photon number resolving detectors. The initial Gaussian 
mixed-entangled states are generated by transmitting a two-mode squeezed state through a lossy 
bosonic channel, corresponding to the primary source of errors in current approaches to optical 
quantum communication. We provide explicit formulas to calculate the entanglement in terms of 
logarithmic negativity before and after distillation, including losses in the channel and the photon 
detection, and show that one-copy distillation is still possible even for losses near the typical fiber 
channel attenuation length. A lower bound for the transmission coefficient of the photon-subtraction 
beamsplitter is derived, representing the minimal value that still allows to enhance the entanglement. 

PACS numbers: 03.67.Mn,03.67.Hk, 42.50.Dv 



I. INTRODUCTION 

Entanglement of composite systems represented by 
continuous quantum variables is of conceptual impor- 
tance for studying fundamental questions of quantum 
mechanics [l[ and promises to be useful for potential real- 
world applications in the fast developing field of quantum 
information processing 043- However, in general, en- 
tanglement is a fragile resource which is easily degraded 
when it interacts with an uncontrollable environment, 
for example, in a communication channel in form of a 
lossy and noisy, optical fiber. In order to overcome this 
problem of entanglement degradation, typically, proto- 
cols such as entanglement purification or distillation will 
be utilized, as originally proposed for qubits IS]. These 
schemes detect errors and are usually probabilistic, as op- 
posed to deterministic approaches such as quantum error 
correction. More generally, entanglement distillation can 
be defined as any scheme that creates one or more en- 
tangled pairs of higher entanglement from one or more 
copies of initially imperfectly entangled pairs by means 
of local operations and classical communication. 

Although many impressive experiments for the distil- 
lation of pure or mixed discrete-variable entanglement 
have been reported [H-IH, distilling continuous- variable 
entanglement appears to be rather different and, in gen- 
eral, harder to achieve. The difficulty arises mainly from 
the necessity of a non-Gaussian element for distilling 
Gaussian entangled states [lOMl 2| | . For instance, in order 
to distill the quantum optical two-mode squeezed state 
(TMSS) whose quadratures obey Gaussian statistics, one 
must introduce at least one non-Gaussian operation, in 
form of a non-Gaussian ancilla or a non-Gaussian mea- 
surement. The so-called photon subtraction (PS) strat- 
egy, first introduced by Opatrny et al. [H[ , is one of the 
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experimentally most readily available operations beyond 
the Gaussian regime. It enables one to modify the Gaus- 
sian statistics of a given TMSS and therefore serves as 
a possible approach to entanglement distillation of such 
Gaussian states The basic principle of the PS tech- 
nique is very simple and can be implemented using a 
beamsplitter and photon measurements. 

After Opatrny et al. 's pioneering work [l3j , many ef- 
forts have been made to improve the performance of 
such an entanglement distillation protocol. Olivares et 
al. [1 51 ] proposed an inconclusive PS method, which em- 
ploys a more realistic on-off photon detector in order to 
enhance the entanglement. Kitagawa et al. [l6j pre- 
sented a fairly complete theoretical analysis of this type 
of distillation, including a numerical evaluation of the 
entanglement before and after distillation. Moreover, a 
multi-mode theory for frequency mode matching in the 
photon-subtracting operation was derived [17| . Even an 
experiment implementing Opatrny's method has already 
been reported [18] , These efforts, both on the theo- 
retical and the experimental side, are examples for the 
more recent attempts to combine discrete-variable and 
continuous-variable approaches to optical quantum in- 
formation processing [19(. 

The original scheme by Opatrny et al. as well as its 
theoretical refinements and extensions all refer to a single 
copy of a pure, Gaussian entangled state which is distilled 
into a non-Gaussian entangled state. This kind of distil- 
lation is sometimes referred to as entanglement concen- 
tration, distinct from entanglement purification proto- 
cols in which initially mixed-entangled states are purified 
and thereby turned into states with higher entanglement. 
Even though usually such entanglement purification is 
applied to two or more copies of entangled states 
one mixed-entangled copy may also be distilled through 
local, generalized measurements, similar to those for con- 
centrating a s ingle p ure-state copy into a maximally en- 
tangled state |20l - [2a |. In the mixed-state case, however, 
both parties sharing the initial state must perform local 
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measurements and communicate their results. 

In this paper, we provide a detailed analysis for the 
one-copy distillation of mixed continuous-variable entan- 
glement, using beamsplitters and experimentally feasible 
photon detection techniques. In other words, similar to 
those one-copy schemes mentioned in the preceding para- 
graph, we shall consider non-Gaussian, generalized mea- 
surements locally performed on the two modes of the 
initial Gaussian state. Note that PS ideally corresponds 
to maps like a\n) = y/n\n — 1) or a\a) = a\a), with the 
photon annihilation operator being non-unitary and the 
resulting states being not normalized - a manifestation 
of the probabilistic nature of the PS process. The corre- 
sponding generalized measurements are realized through 
beam splitter transformations, locally acting on the two 
signal modes and additional ancilla vacuum modes, and 
subsequent photon measurements. Note that various ex- 
periments have already demonstrated how versatile PS 
is for non-Gaussian state engineering [2J-126J. Local fil- 
ters for entanglement concentration of one copy of a pure 
TMSS were previously considered using Kerr interactions 
or CQED [27]. Multi-copy distillations of noisy versions 
of a TMSS have been proposed as well [2lj-[30l|. 

A further recent one-copy scheme for distilling mixed 
continuous variable entanglement was proposed in 
Ref. [3lJ. Different from our analysis, this proposal [3l[ 
employs detections of collective excitations in atomic en- 
sembles for the non-Gaussian operations. Further, it uses 
only operational entanglement measures, namely telepor- 
tation fidelities. Here, we shall calculate both fidelities 
and, in particular, logarithmic negativities for the dis- 
tilled states. 

While in the theoretical analysis of Ref. [l6j], the ini- 
tial states are pure (and become mixed only for the case 
of on-off detections), the input states of the recent ex- 
periment [18| were mixed due to experimental imperfec- 
tions such as the complication to prepare a perfectly pure, 
minimum-uncertainty squeezed state. Although even the 
idealized pure-state versions of that experiment slightly 
differ from the TMSSs used in our analysis (by local 
squeezing operations), the present article also provides 
a more general theoretical foundation of the experiment 
described in Ref. [HI]. It proves the possibility and fea- 
sibility of realistic PS-based distillation of a TMSS, even 
in the presence of high losses. At the same time it illumi- 
nates the applicability of Opatrny's PS distillation pro- 
tocol [13| and provides more details on how to improve 
entanglement in a general, realistic mixed-state scenario; 
for instance, in optical quantum communication using 
lossy fibers of nearly one attenuation length as poten- 
tially used in a quantum repeater (32j |. Moreover, some 
of our results are fully analytical, thus further extending 
the theory presented in Ref. [16]. 

The paper is organized as follows. Firstly, in Sec. II, 
we will give a description of our scheme for entanglement 
distillation, along with the method for generating of our 
mixed entangled state. In Sec. Ill, the definition of log- 
arithmic negativity is briefly summarized. With such a 



figure of merit for entanglement, the amount of entan- 
glement before distillation is explicitly derived. Sec. IV 
is devoted to the entanglement distillation with conven- 
tional on-off detectors. An exact analytical formula for 
the entanglement distillation of pure TMSSs is obtained, 
which previously was only numerically investigated by 
Kitagawa et al. [l6|. In Sec. V, we discuss two different 
strategies of applying photon number resolving detectors 
(pure and mixed PNR detectors). The success probabil- 
ity of distillation and the corresponding lower bounds Tl 
are studied, respectively. In Sec. I VII to further illustrate 
our results, we calculate an operational measure of entan- 
glement (the fidelity in quantum teleportation) , leading 
to yet another way to compare the entangled states be- 
fore and after distillation. Finally, we conclude with a 
discussion and summary in Sec. VII. 



II. LOSSY BOSONIC CHANNEL AND PHOTON 
SUBTRACTION 

Let us first introduce the amplitude-damped TMSS 
with which we are mainly concerned in this paper. This 
kind of mixed-state entanglement can be obtained by 
sending each mode of a pure TMSS through a lossy 
bosonic channel. In our theoretical analysis, we shall sim- 
ulate the lossy channel through an extra vacuum mode 
and a beam splitter. This is the simplest model for 
mimicking realistic fiber-optical light propagation, where 
more and more signal photons are gradually absorbed on 
the way during the channel transmission. 

The entire scheme for our distillation protocol is shown 
in Fig. [TJ The initial pure TMSS is given by 

OO 

\4>)ab = ^2a n \n) A \n) B , a n = V 1 - A 2 A", (1) 

with A = tanh(r) representing the degree of squeezing 
and A, B referring to the two transmitted modes. Two 
beamsplitters with transmission coefficient Tq — r\ and 
auxiliary vacuum modes E, F are put into the ideal chan- 
nels in order to simulate two lossy channels with trans- 
mission efficiency rj. The PS-based distillation is imple- 
mented using two more beam splitters (each with trans- 
mission T) and photon detectors. Due to the finite re- 
flectance of these beam splitter, the photon subtraction 
process heralded by the photon detectors is a probabilis- 
tic process and successful distillation may occur whenever 
both detectors register non-zero counts. 

To give a systematical analysis, in this paper, we 
shall consider entanglement distillation with two different 
types of detectors, namely conventional on-off detectors 
(e.g., avalanche photondiodes, APD) and photon number 
resolving detectors (PNR) , both of which are widely used 
in quantum optical experiments. 
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FIG. 1: Scheme of continuous-variable entanglement distilla- 
tion. The initial state \ip) ab is a pure two-mode squeezed 
vacuum state, both beams of which are transmitted through 
beamsplitters with transmittance To = r\ in order to simulate 
a lossy bosonic channel of transmission rj. The beam split- 
ters with transmittance T including the photon detectors are 
used to achieve the photon subtractions for entanglement dis- 
tillation. The input states of the C, D, E, F-modes are pure 
vacuum states. An event of successful distillation is heralded 
when both detectors register non-zero counts. 



III. LOGARITHMIC NEGATIVITY AND 
ENTANGLEMENT BEFORE DISTILLATION 

Following the definitions of Ref. [l6j], we will use the 
logarithmic negativity as a figure of merit to quantify en- 
tanglement. The logarithmic negativity [33j-|35[ is a rela- 
tively easily computable measure of entanglement; more 
precisely, it is an entanglement monotone, both under 
local operations and classical communication and under 
positive partial transpose preserving operations. 

The logarithmic negativity of a bipartite state p AB = p 
is defined by 



E N (p) = log 2 (l + 2N(p)), 



(2) 



in which N{p) is defined as the negativity given by the 
absolute value of the sum of negative eigenvalues of the 
partially transposed density operator p TA . Here and 
throughout, without loss of generality, we will perform 
the partial transpose operation with respect to the A- 
mode. 

We shall now quantify the amount of entanglement of 
the amplitude-damped TMSS. Note that this state in- 
cluding the damping effect remains a Gaussian state, and 
hence its logarithmic negativity could be directly com- 
puted from its second-moment covariance matrix through 
the corresponding symplectic eigenvalues of the partially 
transposed state 0tJ|- However, for our purposes, since 
the PS process will lead to non-Gaussian states, it is more 
useful to achieve a more general entanglement quantifi- 
cation expressed in the photon number basis. This is 
similar to the approach of Ref. [r| , but with the impor- 
tant distinction that our states are mixed states from the 
beginning, both before and after the distillation. 



First of all, let us denote the beamsplitter coupling 
between modes k and I as [l6j |. 



V k i(0 o ) = exp 



h(a k ai - a k aj) 



(3) 



with 9q = arctan(y / (1 — rj)/rj) and a^m being the pho- 
ton annihilation operator of the fc(7)-mode. The unitary 
state evolution before entanglement distillation can be 
formulated as follows, 



= V AE (0 O ) ® V BF (0oM) AB |0} B |0) F , (4) 



where |0) E , |0) F are the initial vacuum states of the loss 
modes. Through direct calculation, we have 



\^)abef = E E a nUk£,ni\n-k) A \n~l) B \k) E \l) F , 

n=0 k,l=0 



where m = 0, 1, • • • ,n and ( n ) = — rr^ — rr is the bino- 

' ' ' \mJ m!(n— my. 

mial coefficient. 

The mixed state p AB before entanglement distillation 
is obtained by tracing over the loss modes E and F, 

p mlx = p AB = Tr FF [|*) ABFF (*|] 



oo n m 



E EE/-* 

m,n— z—O 7—0 

\n-j) B {m-j\ 



n — i) ,(m — i\ 



(6) 



with fnmij being a real positive coefficient, fnmij = 

(XnQ-m £,ni ^mi £,nj £,mj ■ 

Similar to the case of a pure TMSS [H, [36| , the partial 
transpose of the density matrix ^ is still block diagonal 
in the photon number (Fock) basis. We have 

oo n m 

P T a A b = E EE f n ' m 'ij \ m ~i,n- j) AB (n -i,m-j\ 

m,n=0 2—0 j=0 
oo K 



= E^j^-^bO^-jI. 

K=0 i,j=0 



(7) 



with 



a 



(K) 

hi 



(1-A 2 ) 



n 



K oo 



x (A — At?) 



' / n—no V v / v 

(i+j+2n) (n + iy.(n + jy 



KM{n + i + j - K)V 
n = max{0, K — i — j}. (8) 

Thus, the negativity of p m ix can be equivalently ob- 
tained by solving the eigenvalue problem of a series of 
(K + 1) x (K + 1) sub- matrices 



C 



(K) 



i=0, — ,K, j=0, --,K 



(9) 
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for K = 0, 1, • • • , oo. Such a block sub-matrix method is 
quite efficient in numerical evaluations [1 61 ] . Indeed, the 
matrix Ck has a very useful symmetry property which 
will finally simplify the whole problem. 
Theorem 1. Ck is a double symmetric, i.e., both sym- 
metric and centrosymmetric matrix. 
Proof. Symmetric property follows directly from the i,j 
exchange invariance in the definition Eq. ([8]). Therefore 

we only need to prove — C^} i K _j for arbitrary i, j. 
Now consider any i,j such that for i+ j < K, we always 
have (K - i) + (K - j) > K. Using Eq. ©, it directly 
follows no(i,j) = K — i — j > and no(K — i, K — j) = 0. 
By replacing the index n = n' + no(i, j) in the summation 



of C[j and noticing (^) = (^ij, the relation C 



C 



(K) 



i,3 



K-i,K-j 



can be straightforwardly obtained. 



Theorem 2. The negativity of p m ix can be uniquely 
determined by the skew diagonal entries of matrix Ck: 



N{p mix ) = \{fl ^[JkCk] - 1^ 



(10) 



with Jk being the anti-identity matrix, 



/0 
' 



i,j=0.--,K 



1 
Vl 



• !\ 
•10 



• 

• o/ 



(11) 



Proof: According to Ref. [33], for a (K + 1) x {K + 1) di- 
mensional matrix Ck, there always exists an orthogonal 
matrix J7, such that for K odd, 



UCkU 1 



and for K even, 



A — JM 
A + JM 







(12) 



UC K U 7 



A - JM 

q 
V2x A + JM 



(13) 



where A, M, and J are each L "2^J x L "^^J matrices 
with elements 



A, 



- C [K) 



M,; 



C {K) 



Ji.H = S, 



for i,j = 0,1,---, [(K - 1)/2J, and [---J and [• • ■] are 
the floor and ceiling functions, respectively. 

This shows that the eigenvalues of Ck are the same as 
the eigenvalues of A — JM and A + JM in the case of K 

odd and the same as A — JM and [ % , ^ ) in 

V V2x A + JM J 

the case of K even. 

Moreover, it can be always checked that the sub-block 
A — JM contains all the negative eigenvalues of matrix 
Ck- In fact, for the matrix Ck as defined in Eq.©, 



A— JM is always negative-definite, whereas the sub-block 

A + JM and \ f- are positive-definite. 

V V2x A + JM J P 

Thus, the absolute value of the sum of the negative eigen- 
values of C K is given by \Tr[A - JM]\ = Tr [JM — A] = 
|Tr[JKCK — Ck]- The negativity of the whole matrix 

Pmix follows as 

^ oo 

N( Pmix ) = - ^ Tr[J K C K - C K ] 

K=0 

= ^i>[J K C K ]-^, ( 14 ) 

where in the second line, we have imposed the normal- 
ization condition 



^Tr[C K ] = Tr(^ ) = T r(^J = l. 



(15) 



K=0 



Thus, following the definition in Eq. the logarith- 
mic negativity of the state in Eq. ^ can now be easily 
calculated as 



E N (fimix) = log 2 Tr[J K C K ] 
1 + A 



= log 2 



1 - A(2?7 - 1) 



(16) 



By setting rj = exp(— jt), our result agrees with that pre- 
sented in Ref. [3g]; however, our derivation leads to a sim- 
ple, closed expression as a function of the input squeezing 
and channel loss. For this still Gaussian state before dis- 
tillation, we can also confirm our result by calculating the 
symplectic eigenvalues of the partially transposed state. 

Following our formalism above, we merely need to cal- 
culate the skew diagonal entries {C^_ i }j = o,...,K hi or- 
der to obtain the logarithmic negativity for the two-mode 
mixed entangled state. It is important to note that Theo- 
rem 1 and Theorem 2 can also be applied to calculate the 
logarithmic negativity after entanglement distillation. In 
fact, PS on both transmitted modes does not change 

the symmetry and centrosymmetry of the partially trans- 
r 

posed density matrix p A B, provided that both detectors 
obtain the same measurement results. The only difference 
is that the state after PS is not normalized. One should 
then specify the normalization factor (i.e., the trace of 
Pab = Ejf ^ [^k]) for the different types of detectors 
and the different detection strategies. This enables us 
to extend the analytical formulas for continuous- variable 
entanglement from pure states to mixed states, including 
the Gaussian state p m ix before distillation as well as the 
non-Gaussian states after distillation using on-off detec- 
tors or PNR detectors (in pure and mixed strategies, see 
below) . 
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IV. DISTILLATION USING ON-OFF 
DETECTION 

For convenience, let us first give a general descrip- 
tion of photon detectors. Suppose the detector can re- 
spond with M, different measurement outcomes. Ac- 
cordin g^ to the theory of generalized quantum measure- 
ments [38l l39j , such a measurement device can be com- 
pletely characterized through a set of positive-definite 
operators {IIfc|fc = 1,2, ■• • , M}, corresponding to a pos- 
itive operator- valued measure (POVM). The quantum 
measurement is probabilistic: for a given input state 
g, the probability that the detector gives outcome k is 
Pfe = Trfllfcg]. The condition that the total probabil- 
ity is normalized corresponds to Y2k=i ^fc — 1> with 1 
representing the identity operator. 

The photon detectors usually employed in quantum op- 
tical experiments, such as avalanche photodiodes (APD) 
operating in the Geiger mode, correspond to a measure- 
ment device with only two measurement outcomes: off 
(no photons detected) and on (one or more photons de- 
tected). Expressed in the Fock basis, the positive opera- 
tor description of an ideal on-off photon detector is then 
given by jfiW), f[( on )}, with 



n(°/7) = 

n(°") = i_n(°//) = j2\k){k\. (17) 



Based on the formalism and the notations above, we 
can now proceed with the entanglement distillation pro- 
tocol in Fig|T] Assuming that the two beamsplitters for 
PS have the same transmittance T (reflectance coefficient 
R = 1 — T), the state evolution of the whole PS process 
can be described by 



Pabcd 

p(on, on) 



= V [ Pmix ® |0) c (0| <8 |0) D <0|] Vt, (18) 



,1 



£[(on) (g, n(on) 



P(pn, on) 



where V = V AC {9)®V BD (6>), = y/R/T, and p(on, on) is 
the normalized output state; P(on, on) is the probability 
of detecting non-zero photons in both detectors, 



P(on, on) = Tr 



f[(°n) n(on) 



(19) 



where this time the trace is over all four modes ABCD. 
Using the same method as in Sec. IIII1 analytic formulas 
for the entanglement and the success probability can now 
be derived. The unnormalized, partial transpose p T jf B is 
again block diagonal with respect to the -FT-subspaces. 
We have 



fc=i 



C-j (on, on) 



n = max{0,K + 6 - i - j - 7}. (20) 



The probability of success can be evaluated as where we define T = 1 — r/i? and R=l — rfT. 

00 K 



P(on, on) = C iT ( on ' on ) 

K"=0 i=0 

A 2 (l-f) 2 (l + A 2 T) 
(1-A 2 T)(1-A 2 T 2 ) : 



After state normalization, the logarithmic negativity 
(21) can be also analytically obtained: 



E N (p(on,on)) = log 2 



(1 - X 2 )nR 



log 2 



l-Ar ? r) 2 -A 2 (l-? ? )i? 

(1-A 2 f)(l-A 2 f 2 ) 
(1-T) 2 (1 + A 2 T) 



log 2 



R 



1 — 77 



(1 - A)(l - \{2nT - 1)) (1 - XnT) 2 - A 2 (l - n) 2 



(22) 
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In the following discussions, to be more specific, we 
shall choose two typical values for the channel transmis- 
sion rj in order to study the entanglement properties after 
distillation. 



A. Pure TMSS: n = 1 

In the literature, PS-based distillation of a pure TMSS 
has already been numerically treated in Ref. [1 61 ] . In that 
work, due to the extremely intensive numerical compu- 
tation for diagonalizing a large square matrix, only the 
low-squeezing regime A < 0.9 was investigated and high 
photon number terms were neglected. However, based 
upon our analytical result in Eq. (|22j) . the performance 
of entanglement distillation in the large-squeezing (high 
photon number) regime 0.9 < A < 1.0 can also be exam- 
ined. 

In fact, by simply setting rj = 1, we obtain 



E N (p) v=1 = log. 



(1 + A)(1- A 2 T)(1 + AT) 



P(on, on)^ 



(1-AT)(1 + A 2 T)(1 

A 2 (1-T) 2 (1 + A 2 T) 
(1- A 2 T)(1-A 2 T 2 )' 



AT 



XR) 
(23) 

(24) 



Surprisingly, for a given beamsplitter with finite trans- 
mission coefficient < T < 1, the output entanglement 
exhibits non-monotonic dependence of the initial squeez- 
ing parameter A. The finite transmission coefficient of 
the beam splitter has a degrading effect on the output 
entanglement. When A — > 1, a pure TMSS has infinite 
entanglement. However, when one uses the beamsplitter 
together with on-off detectors to implement the distil- 
lation, one will always get finite entanglement. In fact, 
the optimal squeezing parameter A (referred to as X op t) 
which maximizes E^(p) is strictly smaller than 1. This 
result is certainly of experimental significance in order to 
optimize the distilled entanglement: it may not be nec- 
essary to prepare as much initial squeezing as possible to 
maximize the final entanglement; some finite-squeezing 
value will be optimal. 

In Fig. [5] (a) , we show the logarithmic negativity of the 
distilled TMSS (rj = 1) for different beamsplitter trans- 
missions (T = 0.1, 0.5,- ••, 0.99). In Fig. [2] (b), we 
give a plot of the probability of successful distillation, 
P(on, on). Figured] (c) shows the optimal X op t as a func- 
tion of T, while Fig.[2](d) presents the maximal En(p) at 
X = X op t- Even with infinite squeezing and a non- lossy 
channel, n = 1, we cannot approach infinite entangle- 
ment after distillation. In fact, when A — > 1, in Eq. (I2U1) . 
the logarithmic negativity scales as £jv(/o) = log 2 jzTf ■ 

In the above distillation protocol, there exists a non- 
trivial lower bound Tl for the transmission coefficient T 
below which the PS scheme based on beamsplitters and 
on-off detectors ceases to improve the entanglement. In 
Fig. [2] (a), it is shown that the distillation protocol ef- 
fectively no longer works for T — 0.10, 0.50. The entan- 



glement after distillation is then even smaller than that 
before distillation. Indeed, requiring £jv(/5) > E]y(p m i X ), 
the transmission coefficient T satisfies Tl < T < 1, with 



A(A 2 



A-l) 



2Vesn 



1 

3A3 

A 2 (A 4 + 2A 3 -4A 2 +4A + 1), 
A(A 3 + 8A 2 -3A + 2), 

'3A 3 C-2A(A 2 + A-1)£ 



arccos 



(25) 



The quantity Tl in Eq. (|25|) is a monotonically increas- 
ing function of the squeezing parameter A. When A — > 0, 
we have Tl = 1/2. In the other extreme case, when A 
approaches 1, it follows that Tl — > 1. This, on the other 
hand, proves the degrading effect of the transmission co- 
efficient T: in the high photon number regime (especially, 
for A — > 1), any finite transmission < T < 1 is smaller 
than Tl = 1 and the entanglement of the state after 
distillation, as illustrated in Fig. [TJ is finite and hence 
smaller than the infinite entanglement before distillation. 
We give a detailed description of the behavior of Tl in 
Fig.H] 



B. 3dB amplitude-damped TMSS: n = 0.5 

In Fig. |3] we show the success probability and the 
logarithmic negativity of the amplitude-damped TMSS. 
Here, the amplitude-damping process is simulated by a 
3dB beamsplitter: r\ = 1/2. Similar to the distillation 
of a pure TMSS, the entanglement of the distilled state 
including amplitude damping is again degraded by the 
finite transmission coefficient. There also exists a finite 
X op t with < X op t < 1 which maximizes the output en- 
tanglement. At the same time, the lower bound Tl for 
the transmission coefficient still increases monotonically 
from Tl = 1/2 to Tl = 1, when A varies from to 1. 

In Fig. 21 we show a plot to describe the relation be- 
tween Tl and A, for rj varying from 0.01 to 1. It is shown 
that for larger channel losses (smaller rj), more trans- 
missive beamsplitters (larger T) are needed in order to 
achieve distillation. Furthermore, for all < rj < 1, the 
values Tl vary from 1/2 to 1, which means a beamsplit- 
ter transmission T > 1/2 is a general, necessary condi- 
tion for distilling amplitude- damped TMSSs using on-off 
detectors. However, in Sec. [V] we will find that such a 
necessary condition can be circumvented by employing a 
more sophisticated detection strategy, for instance, using 
photon number resolving detectors. 



V. DISTILLATION USING PHOTON NUMBER 
RESOLVING DETECTION 

In quantum communication and computation, using 
photon number resolving detectors may lead to vari- 
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FIG. 2: (Color online) Performance of distilling a pure TMSS (rj = 1) with beamsplitters and on-off detectors, (a) logarithmic 
negativity of the output state for T = 0.10,0.50,0.70,0.90,0.99, respectively. The dashed line corresponds to the logarithmic 
negativity of the TMSS before distillation (Eg. ll 161) with 77 = 1) . The red squares in each curve indicate the maximum values of 
En{p)- (b) success probability, i.e., the probability that both detectors record the "on" results (Eg. (|24|) '). (c) X op t as a function 
of T (see text for more information), (d) Maximal value of En(p) at A = \ op t- 
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FIG. 3: (Color online) Performance of distilling a 3dB amplitude-damped TMSS (7/ = 1/2) with beamsplitters and on-off 
detectors, (a) logarithmic negativity of the output state for T — 0.10,0.50,0.70,0.90,0.99, respectively. The dashed line 
corresponds to the logarithmic negativity of the amplitude-damped TMSS before distillation fEa.. (|16[0 . The red squares 
indicate the maximum values of En(p)- (b) success probability of distillation for various T (Eq. (|24p ). (c) \ opt as a function of 
T. (d) Maximal value of En(p) at A = X opt . 
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FIG. 4: (Color online) Lower bound Tl for distilling an 
amplitude-damped TMSS using beamsplitters and on-off de- 
tectors. The channel transmissions rj shown, from top to 
bottom, are 0.01,0.1,0.5,1. Tl increases monotonically with 
squeezing A. 



ous important a ppli cations, such as linea r-op tics quan- 
tum computing [40J , quantum repeaters [4l| , quantum 
state discrimination |42j, and quantum superresolution 
[43d . Recently, a photon number resolution of up to 10 
photons was demonstrated 44]. In the following, we 
shall continue investigating PS-based entanglement dis- 
tillation protocols, but we will replace the on-off detec- 
tors by PNR detectors. In our analysis, we will refer to 



two strategies: (1) pure PNR detection strategy and (2) 
mixed PNR detection strategy. 



A. Strategy 1: pure PNR detection 

For simplicity, let us consider a perfect PNR detec- 
tor which has a unique response for every input photon 
number state. The corresponding POVM operator for 
detecting I photons is 



5> 



= 1. 



(26) 



This kind of measurement is pure in the sense that the 
operators IR (t — 0, • • • , 00) are extremal in the convex 
set of all POVMs. Now suppose both PNR detectors in 
Fig. [1] give the same photon number £, then, according 
to Eq. (fl"8"|). the output state can be written as 



p{£,£) 



ABC D AB 



1 



n 



IR 



P(l,t) 



■ (27) 



After direct calculation, we obtain the matrix elements 

r 

of the partially transposed matrix p^ B (unnormalized) 
in the A-subspace, 



rp \K 00 



=(i-A 2 ) 



1-7/ 



j + £\fi + A (K -i + A (K-j + £ 



riQ — max{0, A — i — j}, 



'% +n + A (j + n + A / i + n + £ \ / j + n + l N 
v n J \ n J \i+ j + n — kj\i + j + n — K j 



(28) 



as well as the success probability, 



P{l,i) 



1-A 2 
1 - A 2 T 2 



Xr]R 



1 - X 2 T 2 



21 1 



k=0 



\2k 



(29) 



The logarithmic negativity of the output state then be- 
comes 



E N (p(£,£)) = (2t+l) log 



1 + AT 



1 - X(r]T + 77 - 1) 



loe 



,k=0 



E( l) (A-A^l-Ar/T) 



2i-2k 



log 2 



E 

fc=0 



(AT) 



2 k 



(30) 



In Fig. [5j we show the logarithmic negativity and 
the success probability for distilling a 3dB amplitude- 
damped (77 = 1/2) TMSS. The counted photon numbers 
are t — 1,2,3,4. Compared with the distillation using 
on-off detectors, the PNR-based distillation has the fol- 
lowing characteristics: 

(1) For I > 2, the PNR detectors outperform the on-off 
detectors by a significant amount for small squeezing A. 
The more photons are detected, the higher the entangle- 
ment will be. However, this improvement becomes neg- 
ligible for large squeezing A, for which the lower bound 
T L will be much greater than T = 0.95 (the value used 
in our calculation). 
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FIG. 5: (Color online) Comparison of the performance of en- 
tanglement distillation between PNR detectors (strategy 1) 
and on-off detectors. The green dashed lines correspond to 
the case of on-off detectors. The blue thick dashed lines in 
(a) represent the entanglement before distillation. Logarith- 
mic negativity and success probability are shown for PNR 
detectors with counted photon numbers I — 1,2,3,4. The 
other parameters are set to r\ — 1/2, T = 0.95. The red arrow 
in (b) indicates the regime 0.9950 < A < 0.9999 (3 < r < 5). 
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FIG. 6: (Color online) Lower bound Tl for PNR-based dis- 
tillation (strategy 1) of an amplitude-damped TMSS with 
the numbers of counted photons I — 1,2,3,4. In each plot 
(a)(b)(c)(d), the channel transmittance r\ is chosen to be 
0.01, 0.1, 0.5, 1, from top to bottom, and Tl ranges from 
1/(^+1) to 1. 



(2) The success probability of the pure PNR distilla- 
tion strategy decreases exponentially with the number of 
photons detected in each PNR detector, as can be seen 
in Eq. (j2"9"|) . As a consequence, the probability P(£,£) is 
much smaller than the success probability P(on, on) for 
on-off detectors. To be more specific, we show a plot of 
P(£,£) as a function of r = arctanh(A) in Fig. [5] (b). In 
the high-squeezing regime (0.9950 < A < 0.9999) (3 < 
r < 5), we observe a peak of P(£,£). This is because 
larger squeezing results in more photons in each trans- 
mitted mode (A and B) and therefore leads to more pho- 
tons to be detected by the PNR detectors. However, too 
large squeezing will shift the number of detected photons 
to a much higher level ^ 4, eventually decreasing the 
detection probability for the £ = 1, 2, 3, 4 photon number 
cases. 

(3) The lower bound Tl for the transmittance of the 
beamsplitter is shifted by the PNR detection results. In 
Fig. |51 we show T L as a function of the number of photons 
detected, £ = 1,2,3,4, and the channel efficiency, w = 
0.01,0.1,0.5,1. For £ = 1, the bound T L covers the full 
range between 1/2 and 1, similar to Tl for the on-off 
detection protocol (Fig. HJ). For larger £, e.g. £ — 2,3,4, 
the minimum of Tl (at A = 0) is independent of r\ and 
is shifted to !/(£+ 1), thus circumventing the necessary 
condition T > 1/2 for the on-off detection protocols. 



B. Strategy 2: mixed PNR detection 

To improve the probability of successful distillation, 
we introduce another distillation measurement strategy. 
This time we shall still use photon number discrimination 
with PNR detectors, however, in a mixed PNR strategy. 
Such a strategy is experimentally more feasible than gen- 
eral pure PNR detections and similar experiments have 
already been reported in the context of binary coherent- 
state discrimination [42j . 

To achieve entanglement distillation, we make a post- 
selection of the PNR detection results and define the 
POVM operators 



(31) 



e>m 



Again, successful distillation is heralded when both PNR 
detectors record the "on" result. By taking into account 
the contribution of all multi-photon components £ > m, 
the success probability approaches 1 in the case of infinite 
squeezing (A — > 1). For any m, we have 



p(m) _ ry 



\^Pabcd 



1, 



(i - a 2 ) J2 x2n 



k=0 



■n<™> ®n<™> 



-k) 



(32) 
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FIG. 7: (Color online) Logarithmic Negativity (a) and suc- 
cess probability (b) for the mixed PNR method (strategy 2) 
with m = l,m = 2, m = 3. The lower bound Tl (c) varies 
from l/(m + 1) to 1; the other parameters are chosen as 
■q = 1/2, T = 0.95. 



When m = 1, such a strategy is straightforwardly re- 
duced to the conventional on-off detection method in 
Sec. IIVI However, for large m, the analytic formulas 
for success probability Psucc and logarithmic negativ- 
ity En(J)) become rather complicated and we shall only 
present a numerical comparison for different m values 
in Fig. [JJ We still consider the typical example of 3dB 
transmission r\ = 1/2 and highly transparent beamsplit- 
ters, T = 0.95. As can been seen from Fig. [7J (a), for 
smaller squeezing A < 0.5 (r < 0.5493), a significant in- 
crease of entanglement is obtained. For larger squeezing 
A > 0.5, the mixed PNR detection strategy does not im- 
prove the entanglement very much. The corresponding 



probability P s «£c is shown in Fig. [7J(b). 

Finally, in order to find out for which conditions this 
mixed-PNR protocol can improve entanglement, we also 
systematically vary the T values of the beamsplitters and 
calculate the lower bound Tl (Fig. [7J (c)). Interestingly, 
the Tl bounds are similar to the pure-PNR case. For 
A — s- 0, a transmission of T = l/(m + 1) is sufficient to 
enhance the entanglement. However, as A increases, our 
simulations suggest that a monotonically increasing T is 
required for successful distillation. 



VI. OPERATIONAL MEASURE OF 
ENTANGLEMENT 

In this section, we shall consider quantum teleporta- 
tion of coherent states in order to assess the quality of 
the photon-subtracted entangled states. Quantum tele- 
portation is a protocol in which an arbitrary, unknown 
quantum state can be reliably transferred from a sender 
to a receiver. The crucial resource for quantum telepor- 
tation to outperform classical teleportation is an entan- 
gled state shared by the two participants. Intuitively, 
the more entanglement they share, the higher the tele- 
portation fidelity they can achieve. In other words, the 
teleportation fidelity may serve as an operational mea- 
sure of entanglement [H, [lj| . 

In the following, we consider a teleportation experi- 
ment in which the entangled state after PS-based distilla- 
tion is utilized. We assume that the state to be teleported 
is a coherent state, Oi n = |a)(a|. Standard continuous- 
variable teleportation [45[ consists of three steps: (1) Al- 
ice combines one mode of the entangled state, say the A- 
mode, with the input mode in state <x,- n at a 50 : 50 beam 
splitter; then she measures the quadratures variables 
x- = (xi„, - xa)/V2 and p + = (p in +p A )/v / 2- (2) When 
she obtains the classical measurement results for a;_ and 
she then communicates them to Bob via a classical 
communication channel. (3) Using Alice's measurement 
results, Bob applies the corrsponding displacement oper- 
ation D(—j3) — exp(— j3a) B + f3*a B ),f3 — x_ + ip+ on the 
other entangled mode B. The fidelity between <7j n and 
the final state of mode B is related with the quality of 
the shared entanglement. Unit fidelity requires perfect 
entanglement. 

Mathematically, the joint quadrature measurement on 
the input mode cti n and mode A is equivalent to a het- 
erodyne measurement (acting on mode A), expressible as 

M 



i 



(33) 



where here T denotes the transposition operation. For 
the normalized entangled state p AB , the probability for 
outcome j3 is 



(34) 
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After the displacement operation by Bob, the final state 
in mode B becomes 



1 



Pb = 



D(-P)Tr A {p AB fl A (P) ® 1 b ]DH-0), 



(35) 



with a fidelity given by 

F> =(a\p B \a) (36) 

= W) {a + ^ Tt a[Pab^aW) ® l fl ]|« + 

By averaging over all the possible measurement results 
{3, we obtain the average fidelity 

F = J d 2 f3P(f3)Fp 

= ±J d 2 pTT[p AB D{P)a( n D\p)®D(P)\a){a\D\p)] 

= Tr[^ B P ], (37) 

where we define the bipartite operator Op = 
i / d 2 pD{fi) ® D((3)(af n ® \a){a\)DHfi) ® £>t(/3). Using 
the invariance d 2 B = d 2 (/3 + a), Va, and similar methods 
to those in Ref. [49, we find that 



Or 



oo oo 

EE 

_ff=0 i,j=0 




\i,j}(K-j,K-i\. 



(38) 



Moreover, by noticing that the partially transposed Op 
is block diagonal, we can simplify the fidelity (|37|) as 
follows, 



F = 



^lP r AB O r F ] _Ex= Tr[C K Or(K)] 



Tr(p^) 



(39) 



where Op(K) is the K— sub-block matrix (i|Op(K)|j) = 
(i,K-i\6$\j,K-j). 

Thus, using Ck as defined above, the teleportation 
fidelity can be easily evaluated. For example, for the 
state before entanglement distillation, the matrix Ck is 
given by Eq. ((HJ, and the fidelity becomes 



oo oo 



F, 



^ 

K=0 i,j=0 



(K) 



i\0 F \j) (40) 

(1 + A)(2 - AV + A 2 r? 2 (A + 3) - At?(A + 4)) 
2(2 - 2A?7 - \ 2 r] + \ 2 r] 2 )(l - At?)(1 + A - At?) ' 



Similarly, from the definitions in Eq. (|20p and Eq. 
the teleportation fidelity for the PS-distilled states can 
be obtained, respectively. For example, in comparison 
with the logarithmic negativities calculated in Sec. [V] 
we present a numerical evaluation of the teleportation fi- 
delity for pure PNR-distilled entangled states in Fig. [5] 
For I — 1,2,3, the teleportation fidelity is obviously 
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FIG. 8: (Color online) Fidelity of teleporting an unknown 
coherent state tji n = [«){a| using a pure-PNR-distillcd 
amplitude-damped TMSS. The parameters are r/ = 1/2, T = 
0.95 as in Fig. 



improved in the low-squeezing regime (A < 0.75), in a 
similar way to what we obtained for the logarithmic- 
negativity measured entanglement in Fig. [5] However, 
note that the logarithmic negativity is known to have an 
operational meaning (quantified by the quality of quan- 
tum correlations used in quantum teleportation) only 
for symmetric Gaussian states. Indeed, our amplitude- 
damped TMSSs do belong to the class of symmetric 
Gaussian states. However, for the photon-subtracted, 
non-Gaussian states after distillation, the correspondence 
between logarithmic negativity and coherent-state tele- 
portation fidelity is not obvious; even though it is possi- 
ble to relate the second-moment correlations of photon- 
subtracted states with their logarithmic negativities [47| ■ 



VII. CONCLUSIONS 

In conclusion, we have studied a photon-subtraction- 
based entanglement distillation scheme on a single copy 
of a Gaussian mixed state in form of an amplitude- 
damped two-mode squeezed state (TMSS) using beam- 
splitters and various photon detection strategies. The 
photon measurements included on-off and photon num- 
ber resolving detectors, as well as mixed photon number 
resolving detections where the on-off threshold can be 
varied compared to the conventional on-off measurement 
with zero or non-zero photons detected. Exploiting the 
symmetry and centrosymmetry properties of the partially 
transposed density matrix written in the Fock basis, we 
were able to derive explicit formulas for the entangle- 
ment of the non-Gaussian mixed states after distillation 
in terms of the logarithmic negativity. 

We showed that in order to improve the entanglement 
after the imperfect channel transmission of the TMSS 
subject to photon losses, a constraint represented by a 
lower bound for the beamsplitters (used for photon sub- 
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traction) must be satisfied. Our results extend earlier 
work on continuous-variable distillation from pure en- 
tangled states to the more general case of mixed entan- 
gled states, as one usually encounters in most realistic 
situations such as experimental demonstrations [181 ] and 
optical-fiber-based communications. Most importantly, 
even for channel attenuations as large as 3dB, the photon- 
subtraction-based entanglement distillation scheme still 
works fairly well, provided the input squeezing is chosen 
sufficiently small. 

The applicability of our protocol to actual long- 
distance quantum communication, for instance, by build- 
ing up a quantum repeater (32j, depends on various pa- 
rameters. First of all, note that the success probabilities 
in the present scheme are rather low; i.e., as low as or even 
lower than those of the known discrete- variable repeater 
proposals based upon single-photon detections. More- 
over, our results show that for any (sufficiently small) 
initial squeezing A for which the distilled entanglement 
exceeds the input entanglement, there is always a differ- 
ent, effective squeezing value A c fr > A for which the same 
or even higher entanglement can be distributed through 
the lossy channels without subsequent distillations. This 
suggests that our distillation still mainly functions as an 
entanglement concentration, similar to what can be ob- 
tained for photon-subtraction-based distillation of pure 
states. It is important to see that distillation still works 
for mixed states, however, in a potential application, it 
may still be better to use large squeezing from the be- 
ginning without distillation. In this case, the question 
arises how large this input squeezing must be to beat the 
distillation-based protocol. 



More specifically, using our formulas, one can find that 
the logarithmic negativities before and after distillation 
are related by lim A ^i ^ eforo (?7, A) > £|f cr (?7, A ,T), for 
all initial squeezings Ao, all channel transmissions 77, and 
all photon-subtraction transmittances T < 1. Nonethe- 
less, for example, with 3dB losses in the channel (cor- 
responding to an elementary distance in a quantum re- 
peater of almost one attenuation length), the same en- 
tanglement as for transmitting an almost lOdB-squeezed 
TMSS without distillation can be obtained through 
photon-subtraction-based distillation of a roughly 6dB- 
squeezed TMSS after transmission. However, the former 
approach would be deterministic, whereas the latter is 
highly probabilistic, leading to further complications in 
a full quantum repeater such as the need for sufficient 
quantum memories. Further extensions of our scheme, 
including more general measurements and local opera- 
tions on a single Gaussian mixed state or multi-copy dis- 
tillations may prove superior to the protocol presented 
here. 
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